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Two-body scattering of neutral particles in a trap is studied theoretically. The control of the initial 
state is realized by using optical traps. The collisions inside the trap occur repeatedly; thereby the 
effect of interaction can be accumulated. Two periodic phenomena with a shorter and a much longer 
period, respectively, are found. The latter is sensitive to the interaction. Instead of measuring the 
differential cross section as usually does, the measurement of the longer period and the details of 
the periodic behavior might be a valid source of information on weak interactions among neutral 
particles. 
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The determination of interactions among microscopic 
particles is an important topic in physics. Historically, 
the measurement of the cross sections of 2-body scat- 
tering provides an important source of information. For 
neutral atoms and molecules the determination is more 
difficult because the interaction is in general weak and 
the initial momentum of the incident particle is difficult 
to be precisely controlled. However, if the scattering 
occurs in a trap, new phenomena previously unknown 
might emerge. Due to the recent progress in the tech- 
niques of trapping atoms (molecules) by using optical 
traps [H-Qi trapped scattering might be eventually ex- 
perimentally realized. In a previous paper a model was 
proposed to study the trapped 2-body scattering theo- 
retically [|J. Instead of a single collision, numerous re- 
peated collisions have been found. Thereby the effect of 
interaction can be accumulated and enlarged, and might 
be eventually detected. This favors the determination 
of very weak interaction. The emphasis of that paper is 
placed on the study of the spin-flip phenomenon and the 
determination of the parameter go (the strength of the 
channel with total spin zero) of the 52 Cr atoms. This 
paper is also dedicated to trapped scatterings; however 
the emphasis is placed on the study of a specific period- 
icity emerging from the repeated collisions. It turns out 
that the associated period is sensitive to the interaction. 
Therefore, in addition to the measurement of cross sec- 
tions, the observation of the special periodic phenomenon 
might be also a valid way for the determination of inter- 
action. The model, related theoretical derivation, and 
numerical results are given below. 

We propose a device containing two deep optical traps 
initially. One is close to the origin, while the other one 
is far away. Each trap provides a harmonic potential. 
Thus the total potential was U p (r) = ^Mcdpflr — a| 2 + 
|r — b| 2 ), where a and b are two given vector (norms 
a 3> b), and M. is the mass of the particle involved. 
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Each trap contains a particle in the lowest harmonic 
oscillator (h.o.) state. The two particles are assumed 
to be identical bosons (the generalization to fermions is 
straight forward), the interaction is assumed to be spin- 
independent, and uj p is large enough so that the particles 
are well localized initially and the overlap of their wave 
functions is negligible. Suddenly the two deep traps are 
cancelled. Instead, a broader new trap located at the 
origin U evo i(r) — ^Aiu> 2 r 2 is created, ui < lj p . Since the 
initial state is not an eigenstate of the new Hamiltonian, 
the system begins to evolve. The evolution is affected not 
only by U eV oi{r) but also by the interaction F(|r,— r 3 -|). 
In what follows the details of the evolution is studied, 
two-body collisions occurring repeatedly are found, and 
the effect of interaction is demonstrated. 

Let htu and \Jhj 'Aico be used as units of energy and 
length. The normalized initial state is 



* / = -=(l + P 1 , 2 )(^)texp[-|(|r 1 -a| 2 + |r 2 -b| 2 )] (1) 



where Px,2 implies an interchange of 1 and 2, and rj — 
uj p /uj. We consider the case that a is lying along the 
X— axis, while b along the negative Z— axis. When R = 
(rj + r 2 )/2 and r = r 2 — it for the cm. and relative 
motions, respectively, are introduced, and the h.o. states 
of the new trap are selected as base functions, the initial 
state can be expanded as 
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where (r) is a normalized eigenstates of the Hamilto- 
nian — 5t;V 2 + \ixr 2 with the eigenenergy 2n + 1 + 3/2. L 
and / are coupled to the total orbital angular momentum 
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J and M. 
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where the Clebsch-Gordan coefficients are introduced, 9r 
and 4>R are the spherical polar coordinates, and so on, I 
must be even for boson systems. 

The new Hamiltonian in terms of R and r governing 
the evolution is 



(4) 



Using ifi^/ 2 ^ as base functions, the eigenstates of H r can 
be obtained via a diagonalization and can be expanded 




V> M (r)=5>^Viy 2) (r) 



(5) 



where i denotes the i— th eigenstate of the I— series in 
the order of increasing energy. The associated energy 
is denoted by Eij. Reversely, <p^/ 2 ^ can also be ex- 
panded in terms of ipij. Then, starting from $ /, the 
time-dependent solution of H evo i is 

( T ) = e - iH "° ' T */ 
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where r = cot. In principle, the above solution is exact 
only if the summation covers infinite terms. However, 
when the interaction is not strong, qualitatively accurate 
solutions can be obtained if the number of base functions 
is large enough. This is shown below. 

We define the time-dependent one-body density from 
*(r) as 



p(n,r) = J dr 2 **(r)*(r) 



(7) 



In order to obtain p, the Talmi-Moshinsky coefficients 
relating two sets of coordinates are introduced as 

[^(R) ^{ 2 \v)]j M 

ENLn'LJr (1) / \ (1) / m / n 

nilm2h 

The analytical form of the coefficients can be found in 
0-0] ■ From Eqs. © and ©, making use of the orthonor- 
mality of the base functions, the integration in Eq. ([7]) is 



easy to carry out, and it is straight forward to obtain the 
analytical expression of p(ri,r). Incidentally, since the 
wave function is symmetrized, the behaviors of the two 
particles are exactly the same. The observation of only 
one particle is sufficient. 

To obtain numerical results as examples, it is first as- 
sumed that a = 2, b = 0, r\ = 2, and the interaction 
contains a stronger repulsive core and a weaker attrac- 
tive tail as 



V(r) 



V , 

-C 6 /r e , 



if r < 0.2 
else 



(9) 



where Vb and Cq are positive numbers, and Vq >• Cq. We 
define K = 2(N+n)+L+l to control the dimension of the 
base. Mostly, the base functions with K < K max = 20 
are adopted in the following calculation. 




FIG. 1: (Color online) p(ri,r) plotted on the X — Z plane. 
t = 0, 7r/6, 7t/3, 7t/2, 57r/6 and 7r, respectively, for (a) to (f). 
The parameters are rj — 2, a = 2, b = 0, Vo — 10, ro = 0.2 and 
ce = — 0.3ro. The units of energy and length in this paper are 
ftw and y/h/Mio, where oj has not yet been specified. Every 
maximum in the panels is marked by a cross. The values 
associated with the inmost contours of (a) to (f) are 0.224, 
0.098, 0.073, 0.063, 0.098 and 0.224. Thus the peaks in (a) 
and (f) are much higher. Therefore, the particles are better 
localized when they separate from each other. The values of 
the outmost contours are from 0.01 to 0.03. 

When Vb = 10, Cq = —0.3, and r is given at a number 
of values in the early stage of evolution, p(ri,r) plot- 
ted on the X — Z plane is shown in Fig. [TJ [TJl is for 
the initial case. When the evolution begins, the outside 
particle moves toward the center and collides with the 
target particle as shown in [T]d and [T};. When r = ir/2 
both particles are close to the center as in[TJi. Afterward 
a particle begins to leave as in QJ, and will arrive at the 
opposite end at r = n as in [If. Then the process re- 
peats but in reverse direction. When r = 27r the system 
recovers its initial status. This is a kind of periodic be- 
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havior originating from the harmonic trap. If the interac- 
tion is neglected, the recovery would be exact. However, 
due to the interaction, p(ri , r) is not exactly equal to 
p(ri, r + 2ir). With the above parameters, the deviation 
is very small. Say, when ri = a, we have p(a, 0) = 0.2553 
and /?(a, 2ir) — 0.2544. Obviously, in the period (0, 2tt) 
the system undergoes a pair of collisions. When the time 
goes on, a series of head-on collisions occur repeatedly in 
the trap. Although the effect of interaction on a round 
of collision is weak, the effect of many rounds might be 
accumulated and therefore might become strong. This is 
shown below. 



p(a.-0 




FIG. 2: (Color online) p(a, r) plotted against r, where a de- 
notes the initial position of the incident particle and r is given 
only in discrete values 2kn, where k is an integer from are 
60. Three choices of Vo, namely, 20, 10, and 2 are adopted, 
and the associated p are marked with circles, diamonds, and 
triangles, respectively (however the values of p associated with 
an odd k have been neglected just for simplicity). The other 
parameters are the same as in Fig. [1] 

In Fig. [2] p(a, r) is given at r = 2kn, where k is an 
integer. It implies that p is observed at the initial posi- 
tion of the incident particle repeatedly. If the interaction 
is removed, all the symbols in the figure would lie along 
a horizontal line (implying an exact periodicity with a 
period 2-7r). However, the interaction causes a deviation. 
The deviation would become larger if the interaction is 
stronger (the black circles to be compared with the trian- 
gles) and/or if r is larger. It is found that, when the time 
goes on, p(a, 2kir) against increasing k will first arrive at 
a minimum, then arrive at the second maximum which is 
a little lower than the first maximum at r = 0, then again 
a minimum, and afterward arrive at the third maximum 
with a height close to the first maximum. This behavior 
will repeat again and again. E.g., for Vo — 20 (black cir- 
cles), the first, second and third maxima appear at r = 0, 
54-7T, and 1107T, respectively. Whereas for Vq = 10 (red 
diamonds), they appear at 0, 1047T and 2107T. 

On the other hand, making use of the expansion 
Eq. ©, we calculate the overlap | <\I> (0) | >3> (2/ctt)) | which 
varies with k. When k leads to a minimum (maximum) of 
p(a, 2kir), the overlap is small (close to one). For exam- 
ples, with the parameters for Fig.Q] when t = 487r, 104-7T 
and 2107T associated with the first minimum, the sec- 
ond maximum, and the third maximum, respectively, we 
have |(*(0)|*(t))| = 0.103, 0.918 and 0.994. It is further 



noted that the imaginary part of (^(0) 1^(210^)) is very 
small. Thus, , 5(2107r) is extremely close to the initial 
state. Therefore, from the time-dependent Schrodinger 
equation, we know that what happens during the inter- 
val (0, 2107r) will nearly exactly repeat again in the next 
interval (2107T, 4207r), and so on. It implies the existence 
of another nearly periodic behavior. Thus, there are two 
distinct periodic behaviors. One has a period 2-7T, and the 
other one has a much longer period (say, for the above 
case, the period is 2107r). 
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FIG. 3: (Color online) p(ri, r) plotted along the X— axis. The 
curves "1" to "5" in [3^ have t from to it with a step 7r/4. 
Those of [3b have r from 48tt to 49ir with the same step. The 
parameters are the same as in Fig. \T\ 
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FIG. 4: (Color online) The same as Fig. [3] but the domain of 
r is (104tt, 105tt) ingja, and (210tt, 211tt) ingb. 
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The periodic behaviors can be shown in more detail 
by observing directly the densities. Let r x — xi, where 
i is a unit vector along the X— axis. The distribution 
of p(xi, t) along the X— axis is plotted in Figs. [3J and |H 
where r is given at a number of values. Fig. [3^ describes 
the evolution in the interval (0, tt), where the curve "1" 
is for the initial state. From "1" to "5" r goes from 
to 7r with a step 7r/4. One can see how the two particles 
undergo a round of collision. In fact, Fig. Q] and Fig. [3^ 
describe the same thing except that p is plotted on the 
X — Z plane in the former but only along the X— axis 
in the latter. When r goes from 7r to 27r, the process 
occurring in (0, ir) will repeat again but in reverse direc- 
tion. Thereby the cycle with the 2ir period is completed. 
Fig. [3Jd describes the evolution in the interval (487r, 497r) 
associated with the first minimum of the diamonds in 
Fig.[2j The peaks in [3]d are much lower than those of [3ji, 
and these peaks are located at different places. Therefore 
the system behaves differently in the two intervals, and 
the previous clear picture of a head-on collision becomes 
ambiguous. Fig. 0^ is associated with the second max- 
imum. The collision can be roughly seen but is not as 
clear as in Fig. [3^,. Fig. HJd is associated with the third 
maximum, and is nearly identical to |3ji. Therefore the 
cycle with the longer period is completed. 
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FIG. 5: (Color online) Similar to Fig. [2] but Vo is fixed at 
10 and C6 has three choices: — 0.3t"o, — 1.5?"o and — 3rjj. The 
values of p(a, r) are, respectively, marked by black circles, red 
diamonds, and triangles. The other parameters are the same 
as in Fig. [T] This figure has a longer range than that in Fig. [2] 
and only the values of p with r = 8kn are shown. 

In Fig. [51 the strength of the attractive tail has been 
given at three values. In each case a slightly lower peak 
followed by a higher peak (the third maximum) appears 
again. With the three choices of cq, the third maximum 
appears at r = 2107T, 2367T, and 2827T, respectively. At 
these three instants, |(*(0)|*(r))| are all equal to 0.994 
and the associated imaginary parts are very small. It 
implies a nearly exact recovery. Thus the nearly peri- 
odic behavior with the much longer period appears again. 
When 6^0 (i.e., the target particle is not at the center 



initially), the above qualitative features remain. In par- 
ticular, for a specific interaction, the longer period is not 
changed with b. Thus we conclude that the periodic phe- 
nomenon is common to trapping 2-body scatterings. It is 
emphasized that this phenomenon is sensitive to the in- 
teractions. A stronger repulsive (attractive) force would 
lead to a shorter (longer) long-period. 

The accuracy of the above numerical results depends 
on isTmax- As an example selected values of p(a, r) are 
listed in Tab. Uto show the dependence. 



TABLE I: p(a, r) with three choices of K m: 
involved are the same as those for Fig. [T] 



The parameters 



'max 


12 


16 


20 


p(a,0) 


0.253 


0.255 


0.255 


p(a, 50-7I-) 


0.067 


0.068 


0.068 


p(a, 100tt) 


0.229 


0.230 


0.229 



The convergency appears to be satisfying. Thus the 
numerical results obtained by using -Km ax = 20 are accu- 
rate enough in qualitative sense. 

In this paper the traditional 2-body scattering is con- 
sidered under a new environment, namely, in a trap. Due 
to the trap, the two particles collide with each other re- 
peatedly. Therefore, even the interaction is weak; the 
effect of interaction can be accumulated via the repeated 
collisions. Besides, comparing with the case of charged 
particles, the initial scattering states of neutral particles 
are more difficult to control. This disadvantage can be 
overcome by using optical traps. Furthermore, two peri- 
odic phenomena are found. They are essentially caused 
by the trap and by the interaction, respectively. The pe- 
riod of the latter is much longer and is sensitive to the in- 
teraction. Instead of measuring the differential cross sec- 
tion as usually does, the observation of the longer period 
and the details of the time-dependent density might be a 
valid source of information on weak interactions among 
neutral particles. 

The above approach can be easily generalized to the 
cases with various interactions, and/or to the case of 
Fermion systems. It is reasonable to expect that the 
above trapped scattering could be experimentally real- 
ized via the progress of techniques in trapping neutral 
particles by optical traps. 
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